
import io
from os.path import exists
from urllib.request import urlretrieve
from zipfile import ZipFile
from datetime import timedelta
import pandas as pd
import geopandas as gpd

import warnings
warnings.filterwarnings('ignore')


def download_from_zenodo(zip_file, url):
    if not exists(zip_file):
        urlretrieve(url, zip_file)


ais_zip_file = 'data/P1_AIS_Data.zip'
download_from_zenodo(ais_zip_file, 'https://zenodo.org/record/1167595/files/%5BP1%5D%20AIS%20Data.zip?download=1')


NROWS = 10000

def csv_to_df(csv_file_name="nari_dynamic.csv"):
    with ZipFile(ais_zip_file) as z:
        with z.open(csv_file_name) as f:
            df = pd.read_csv(f, nrows=NROWS)
            try:
                df['t'] = pd.to_datetime(df['ts'],unit='s')
                df.drop(columns=['ts'])
            except:
                df['t'] = pd.to_datetime(df['t'],unit='s')
            df = df[['t', 'sourcemmsi', 'lon', 'lat']]
            return df
            
df = csv_to_df()

df.head()


gdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df.lon, df.lat))
gdf


import movingpandas as mpd
mpd.__version__

INFO: Missing optional dependencies. To use the trajectory smoother classes please install Stone Soup (see https://stonesoup.readthedocs.io/en/latest/#installation).

'0.9.rc3'


tc = mpd.TrajectoryCollection(df, 'sourcemmsi', x='lon', y='lat', t='t', min_length=1000)
tc = mpd.MinTimeDeltaGeneralizer(tc).generalize(tolerance=timedelta(minutes=1))
hvplot_defaults = {'tiles':'OSM', 'frame_height':400, 'frame_width':500, 'line_width':5.0}
basic_plot = (tc.hvplot(hover_cols=['sourcemmsi'], **hvplot_defaults) * 
              tc.get_start_locations().hvplot(geo=True, c='white') * 
              tc.get_end_locations().hvplot(geo=True, c='black'))
basic_plot


speed_plot = tc.hvplot(title=f'Speed (m/s)', c='speed', cmap='Viridis_r', colorbar=True, **hvplot_defaults)
speed_plot


(basic_plot + speed_plot).cols(1)


import skmob
skmob.__version__

'1.2.3'


tdf = skmob.TrajDataFrame(df, latitude='lat', longitude='lon', datetime='t', user_id='sourcemmsi')
tdf.plot_trajectory(zoom=9, weight=3, opacity=0.9, tiles='OpenStreetMap', max_users=100)


import trackintel as ti
ti.__version__

'1.1.5'


pfs = ti.io.from_geopandas.read_positionfixes_gpd(gdf, tracked_at='t', user_id='sourcemmsi')
pfs, sp = pfs.as_positionfixes.generate_staypoints(method='sliding')  # dist_threshold=10)
pfs, tpls = pfs.as_positionfixes.generate_triplegs(sp, method='between_staypoints')
tpls.as_triplegs.plot()


import traja
traja.__version__

'0.2.13'


trj = traja.TrajaDataFrame(df.rename(columns={"lon": "x", "lat": "y", "t": "time"}))
trj.traja.plot()

<matplotlib.collections.PathCollection at 0x21c71688280>


trjs = traja.TrajaCollection(df.rename(columns={"lon": "x", "lat": "y", "t": "time"}), id_col="sourcemmsi")
lines = trjs.plot()


help(trjs.plot)

Help on method plot in module traja.frame:

plot(colors=None, **kwargs) method of traja.frame.TrajaCollection instance
    Plot collection of trajectories with colors assigned to each id.
    
    >>> trjs = {ind: traja.generate(seed=ind) for ind in range(3)} # doctest: +SKIP
    >>> coll = traja.TrajaCollection(trjs) # doctest: +SKIP
    >>> coll.plot() # doctest: +SKIP


import transbigdata as tbd


traj_gdf = tbd.points_to_traj(df, col=['lon','lat','sourcemmsi'])
traj_gdf.plot(column='sourcemmsi', categories=df.sourcemmsi.unique().tolist())

<AxesSubplot:>


tbd.visualization_trip(df, col=['lon','lat','sourcemmsi','t'])

Processing trajectory data...
Generate visualization...
User Guide: https://docs.kepler.gl/docs/keplergl-jupyter

KeplerGl(config={'version': 'v1', 'config': {'visState': {'filters': [], 'layers': [{'id': 'hizm36i', 'type': …


from ptrail.core.TrajectoryDF import PTRAILDataFrame
from ptrail.visualization.TrajPlotter import TrajectoryPlotter
import ptrail
ptrail.__version__

'0.7.1 Beta'


pdf = PTRAILDataFrame(data_set=df, latitude='lat', longitude='lon', datetime='t', traj_id='sourcemmsi', rest_of_columns=[])
TrajectoryPlotter.show_trajectories(pdf)

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
Input In [22], in <cell line: 2>()
      1 pdf = PTRAILDataFrame(data_set=df, latitude='lat', longitude='lon', datetime='t', traj_id='sourcemmsi', rest_of_columns=[])
----> 2 TrajectoryPlotter.show_trajectories(pdf)

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\ptrail\visualization\TrajPlotter.py:244, in TrajectoryPlotter.show_trajectories(dataset, weight, opacity)
    241 TrajectoryPlotter._animal = TrajectoryPlotter._create_radio()
    243 # Create the multi selection button.
--> 244 TrajectoryPlotter._selector = TrajectoryPlotter._create_multi_select(TrajectoryPlotter._dataset,
    245                                                                      TrajectoryPlotter._animal.value)
    247 # Create the widgets.
    248 ie = widgets.interactive_output(TrajectoryPlotter._plot, {'value': TrajectoryPlotter._selector})

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\ptrail\visualization\TrajPlotter.py:50, in TrajectoryPlotter._create_multi_select(dataset, animal)
     33 @staticmethod
     34 def _create_multi_select(dataset, animal):
     35     """
     36         Create the multiple selection widget.
     37         Parameters
   (...)
     48 
     49     """
---> 50     dataset = dataset.reset_index()
     52     # Select the animal based on the parameter passed.
     53     to_select = None

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pandas\util\_decorators.py:311, in deprecate_nonkeyword_arguments.<locals>.decorate.<locals>.wrapper(*args, **kwargs)
    305 if len(args) > num_allow_args:
    306     warnings.warn(
    307         msg.format(arguments=arguments),
    308         FutureWarning,
    309         stacklevel=stacklevel,
    310     )
--> 311 return func(*args, **kwargs)

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pandas\core\frame.py:5787, in DataFrame.reset_index(self, level, drop, inplace, col_level, col_fill)
   5785     new_obj = self
   5786 else:
-> 5787     new_obj = self.copy()
   5789 new_index = default_index(len(new_obj))
   5790 if level is not None:

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pandas\core\generic.py:6032, in NDFrame.copy(self, deep)
   5926 @final
   5927 def copy(self: NDFrameT, deep: bool_t = True) -> NDFrameT:
   5928     """
   5929     Make a copy of this object's indices and data.
   5930 
   (...)
   6030     dtype: object
   6031     """
-> 6032     data = self._mgr.copy(deep=deep)
   6033     self._clear_item_cache()
   6034     return self._constructor(data).__finalize__(self, method="copy")

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pandas\core\generic.py:5575, in NDFrame.__getattr__(self, name)
   5568 if (
   5569     name not in self._internal_names_set
   5570     and name not in self._metadata
   5571     and name not in self._accessors
   5572     and self._info_axis._can_hold_identifiers_and_holds_name(name)
   5573 ):
   5574     return self[name]
-> 5575 return object.__getattribute__(self, name)

AttributeError: 'PTRAILDataFrame' object has no attribute '_mgr'


pdf.head(5)


from ptrail.core.Datasets import Datasets
starkey = Datasets.load_starkey()
starkey.head(5)


from ptrail.visualization.TrajPlotter import TrajectoryPlotter as tpl
tpl.show_trajectories(starkey, opacity=0.75, weight=2)


import pymove as pm
import folium as f
import matplotlib as mpl
from pymove.utils import visual
pm.__version__

---------------------------------------------------------------------------
ImportError                               Traceback (most recent call last)
Input In [23], in <cell line: 1>()
----> 1 import pymove as pm
      2 import folium as f
      3 import matplotlib as mpl

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pymove\__init__.py:13, in <module>
     11 from .core.dataframe import MoveDataFrame
     12 from .core.grid import Grid
---> 13 from .core.pandas import PandasMoveDataFrame
     14 from .core.pandas_discrete import PandasDiscreteMoveDataFrame
     15 from .models.pattern_mining import clustering

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pymove\core\pandas.py:47, in <module>
     45 from pymove.utils.log import logger, progress_bar
     46 from pymove.utils.mem import begin_operation, end_operation
---> 47 from pymove.utils.trajectories import shift
     49 if TYPE_CHECKING:
     50     from pymove.core.dask import DaskMoveDataFrame

File H:\miniconda3\envs\geospatial-plus\lib\site-packages\pymove\utils\trajectories.py:23, in <module>
     21 from pandas import DataFrame, Series
     22 from pandas import read_csv as _read_csv
---> 23 from pandas._typing import FilePathOrBuffer
     25 from pymove.core.dataframe import MoveDataFrame
     26 from pymove.utils.constants import DATETIME, LATITUDE, LONGITUDE, TRAJ_ID, TYPE_PANDAS

ImportError: cannot import name 'FilePathOrBuffer' from 'pandas._typing' (H:\miniconda3\envs\geospatial-plus\lib\site-packages\pandas\_typing.py)


def csv_to_pymove_df(csv_file_name):
    with ZipFile(ais_zip_file) as z:
        with z.open(csv_file_name) as f:
            df = pm.read_csv(f)
            return df.head(NROWS)
        
move_df = csv_to_pymove_df("nari_dynamic.csv")

	t	sourcemmsi	lon	lat
0	2015-09-30 22:00:02	245257000	-4.465718	48.382490
1	2015-09-30 22:00:03	227705102	-4.496571	48.382420
2	2015-09-30 22:00:04	228131600	-4.644325	48.092247
3	2015-09-30 22:00:05	228051000	-4.485108	48.381320
4	2015-09-30 22:00:06	227574020	-4.495441	48.383660

	t	sourcemmsi	lon	lat	geometry
0	2015-09-30 22:00:02	245257000	-4.465718	48.382490	POINT (-4.46572 48.38249)
1	2015-09-30 22:00:03	227705102	-4.496571	48.382420	POINT (-4.49657 48.38242)
2	2015-09-30 22:00:04	228131600	-4.644325	48.092247	POINT (-4.64433 48.09225)
3	2015-09-30 22:00:05	228051000	-4.485108	48.381320	POINT (-4.48511 48.38132)
4	2015-09-30 22:00:06	227574020	-4.495441	48.383660	POINT (-4.49544 48.38366)
...	...	...	...	...	...
9995	2015-10-01 02:03:02	227705102	-4.496548	48.382400	POINT (-4.49655 48.38240)
9996	2015-10-01 02:03:03	228051000	-4.485067	48.381332	POINT (-4.48507 48.38133)
9997	2015-10-01 02:03:05	228374000	-4.769649	48.058800	POINT (-4.76965 48.05880)
9998	2015-10-01 02:03:06	227574020	-4.495477	48.383625	POINT (-4.49548 48.38363)
9999	2015-10-01 02:03:07	228138000	-4.775683	48.019770	POINT (-4.77568 48.01977)

Trajectory Visualization Survey¶

MovingPandas¶

Scikit-mobility¶

Trackintel¶

Traja¶

TransBigData¶

PTRAIL¶

PyMove¶

Appendix¶